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Abstract 

Quantum algorithms speeding up classical counterparts are proposed for the problems: 

1. Recognition of eigenvalues with fixed precision. Given a quantum circuit generating 
unitary mapping U and a complex number the problem is to determine is it an eigenvalue 
of U or not. 

2. Given a molecular structure find thermodynamic functions like partitioning func- 
tion, entropy, etc. for a gas consisting of such molecules. 

3. Recognition of molecular structures. Find a molecular structure given its spectrum. 

4. Recognition of electronic devices. Given an electronic device that can be used only 
as a black box how to recognize its internal construction? 

We consider mainly structures generating sparse spectrums. These algorithms require 
the time from about square root to logarithm of the time of classical analogs and for 
the first three problems give exponential memory saving. Say, the time required for 
distinguishing two devices with the same given spectrum is about seventh root of the 
time of direct classical method, for the recognition of eigenvalue - about sixth root. Thus 
microscopic quantum devices can recognize molecular structures and physical properties 
of environment faster than big classical computers. 
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1 Electronic devices and quantum computations 

1.1 Statements of problems and outline of the work 

The aim of this paper is to build effective quantum algorithms for the problems of the following 
types: 

• Given a quantum gate array generating unitary operator U and a complex number uj how 
to determine is it an eigenvalue of U or not (precision of determining eigenvalues is fixed) ? 

• How to recognize a structure of unknown electronic or molecular device given only access 
to its function? 

Here the first problem will be an important intermediate step in the solution of the second 
Q. Consider them sequentially. 

Recognition of eigenvalues. This problem is closely connected with finding of eigenvalues 
distribution or density of states (DOS) that is energy levels Eq < Ei < . . . and dimensions of 
the corresponding subspaces do, di, . . .. DOS plays a key role in calculation of thermodynamic 
functions given by 

F = J2aij)d,e-^ (1) 
j 

for some values a{j) so that the summands quickly converge to zero. Say, if all a{j) = 1 
this expression gives the partition function Q, if a{j) = Ej/Q it gives an average energy, if 

straightforward calculation shows that the simulation of evolution generated by a given Hamiltonian up 
to a time instant r with a fixed accuracy requires of order steps on a quantum computer. This means that 
all results of the paper can be generalized to arbitrary quantum systems. This subject will be elaborated in 
more details in another paper. 
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a{j) = — ^ \n(e~^^^''^^ /Q) - an entropy. Having an efficient method of finding dj we would be 
able to obtain thermodynamic functions and to determine important properties of environment 
consisting of such molecules like heat capacity. The best known classical method of finding 



DOS was proposed by Hams and Raedt in the work |PR|| . Their method requires the time of 
order dimension of the space of states and the memory of the same order (whereas the direct 
method of calculation eigenvalues requires the time of order N^). The ffist quantum algorithm 



for this problem proposed by Abrams and Lloyd in ||AL|| requires the same time 0{N) and 



logarithmic memory. The method proposed in the present work requires the time of order 
square root of classical and memory of order In^ A^. 

The idea of our approach is the following. We shall use combination of Grover search 
algorithm (GSA), revealing eigenvalues by Abrams and Lloyd method ( ||AL|| ) and universal 
quantum function of application App. Abrams and Lloyd method of revealing eigenvalues is 
based on the application of U controlled by ancillary qubit a as 



U, 



cond 




1, 

0. 



Note that it is the direct generalization of Shor's trick which can be obtained if f/ is a multi- 
plication by a given integer modulo q ( |^h|| ). 
Recognition of devices structure. 

We shall tell apart two versions of this general problem: the recognition of molecular struc- 
tures and the recognition of electronic circuits. 

If we want to determine a molecular structure then it is natural to assume that its function- 
ality is given as the spectrum of its Hamiltonian, e.g. the set of its energy levels. Thus here it 
is required to find a quantum system whose Hamiltonian has a given spectrum. 

The problem of recognition of electronic circuits is stated otherwise. An electronic device 
is thought of as a source of electromagnetic fields which can control some quantum system Q. 
Let such a field induces evolution of the system with Hamiltonian H in time frame 6t. Thus we 
have a correspondence: electronic device — > Hamiltonian, 6t. Evolution of a quantum system 
Q induced by this Hamiltonian can be presented as a unitary transform U = e~h^^^. Then, 
given a device C and a value of time t we can associate with it some unitary transformation Uc- 
Assume that we have recognized a circuit C if we find some circuit Ci such that Uc = Uc^ with 
high accuracy. We shall write U instead of Uc for the circuit C that we want to recognize. But 
in fact we shall solve the more general problem when a tested device C can be used as a black 
box acting on n qubits as a function Uc so that if x is an input then Uc\x) is a result of its 
action on this input. Here a tested device can contain its own quantum memory and it can be 
entangled with Q in course of fulffiling U but this entanglement must be then eliminated. The 
existence of such entanglement means that this case cannot be described by the Hamiltonian 
of system Q. For the simplicity we assume that an unknown circuit is built of elementary 
functional elements from some fixed set {^'i, £^2, • • • , E^}. The next natural assumption is that 
a size of circuit is limited by some constant c so that our circuit is some unknown combination of 
c functional elements. Denote by 8 all circuits of the length c. We can encode such C G £^ by a 
string [C] of ones and zeroes so that decoding procedure is easy as well and we can immediately 
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recreate a circuit given its code. Thus we can look through all circuits looking through its 
codes. The same coding can be built for electronic devices. 

A straightforward solution of the problems is clear. For the problem of recognition of 
molecular structures all that we need is to be able to recognize eigenvalues of transformation 
generated by a given circuit. Each eigenvalue of unitary operator has the form ^'^'^ where lo 
is a real number from [0, 1) called frequency. In what follows by spectrum we mean a set of all 
frequencies. Let all frequencies are grouped near points of the form where M is not very 
big, Z = 0,1,...,M — 1. Assume that the acceptable precision of recognition of frequencies is 
1/M. Then having an algorithm for eigenvalue recognition we can apply it again and again 
constructing spectrums generated by all possible circuits and thus find a wanted circuit with 
given spectrum. If we need to recognize a circuit of electronic device we can examine all possible 
circuits taken in some order. Examination of one circuit means that we run it on all possible 
inputs one after another and compare the results with the corresponding result of a tested 
device action. 

For the problem of recognition of molecular structures our method requires the time of order 
sixth root of the time of direct classical method whereas memory saving is exponential. For the 
problem of recognition of electronic circuits our method gives at least square root time saving in 
the case when classical counterparts exist (this is the narrow formulation when a tested device 
generates classical mapping). But in general case an advantage may be more. For example we 
can tell apart two devices with the same spectrum in the time about seventh root of the time 
of naive bruit force. 

To recognize devices quantumly we must be able to store and fulfill operations on the codes 
of different circuits. This possibility is based on the existence of a quantum analog of the 
universal Klini function. This is a unitary operator App such that for all quantum devices C 
and all inputs x App|,T, [C]) = \Uc x, [C]). We assume that for the wide range of quantum 
devices C with c particles C may be encoded as integer [C] in time 0(c) so that the quantum 
complexity of App is 0(c) as well. 

We shall consider here a particular case of the problem when all eigenvalues of U are known 
a-priori or can be obtained beforehand. This restriction is not yet very constraining. To 
illustrate what kind of tasks we shall be able to solve by the proposed method consider a few 
examples of the problem of recognition of an electronic device whose spectrum is known. 

Recognition of quantum algorithms designed as subroutines. Such algorithms must restore 
an input if we apply it twice. Computing a function / they act as \x, b) — > |x', b + f{x) mod 2). 
All known quantum algorithms can be presented in such form. For such quantum algorithms 
their unitary transformation U has only two eigenvalues: 1 and —1. Given a controlling device 
for such algorithm (it may include classical elements and ancillary qubits as well) we can quickly 
recognize its construction. Alternatively, we can quickly find quantum or classical algorithm 
for a given task. 

Consider a "classical" particular case of the recognition problem when U maps each basic 
state to a basic state which means that the matrix of U consists of ones and zeroes and in 
addition U equals U~^. Here evident strategy of recognition takes of order card(£^)2"' steps. 
This case of the problem may be reformulated as the finding of such t that for all s some 
given predicate A{t, s) is true. This is the problem of verification of logical formulas. Quantum 
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solution of it in a time about square root of classical time based plainly on Grover's trick 
was proposed in PCW|| . This method doesn't work in the general case where Uc is arbitrary 



involutive unitary transform, e.g. such that U = U^^. Just this general case is the subject of 
this work. Here we cannot recognize a circuit so easy as in "classical" case because it is difficult 
to compare two quantum states Uc\x) and U\x). 

The general idea of our approach to the recognition of arbitrary electronic devices is the 
following. We shall include a device C which structure we want to recognize into a classical 
controlling part of quantum computer. Thus a tested device generates unitary transformation 
on n qubit system. Then reveal eigenvectors of U using Ucond by the method mentioned above 
and compare them with eigenvectors of circuits from S choosing a circuit giving the best 
approximation. Here GSA will be used in the last step and in the several intermediate steps. 
Assumption about sparse spectrum 

In this paper we shall consider mainly circuits generating sparse spectrums. It means that 
spectrums of operators Uc are so designed that the frequencies are grouped into groups such 
that a minimal distance between frequencies from the different groups is more than 1/M and a 
maximal distance between frequencies from the same group is less than 1/L. For the problems 
of eigenvalues and molecular structures recognition we require that L = IQM that is not yet 
very restricting. For the recognition of electronic devices we shall suppose L ^ M that is 
more strong limitation. Spectrums are called sparse if M = const when — > oo. For sparse 
spectrums our algorithms show the best performance. 

Spectrums that are not sparse are called dense. For dense spectrums our methods give the 
less advantage over classical algorithms (look at section 3.6). An example of dense spectrum: 
ujk = jf, /c = 0, 1,...,A^ — 1. The similar problems for dense spectrums will be studied in one 
of the following papers. 

We write uj' ^ uj iS uj' and uj belong to the same group. For the simplicity assume also that 
for each group of frequencies there exists a number of the form l/M disposed between some 
two frequencies of this group where I is an integer less than M. 



1.2 An abstract model of QC. "Plug and play" technology 

To build algorithms recognizing circuits we need an abstract model of quantum computer (QC). 
QC consists of two parts: quantum and classical. Classical part exactly determines what unitary 
transformation must be fulfilled at each time instant with quantum part and thus plays a role 
of controller for it. These unitary transformations are of two sorts: working transformations - 
which our computer performs itself and query transformations - induced by a tested device: U 

or f/cond- 

We can suppose that a quantum part Q consists of nuclear spins or interacting dipoles (or 
some other quantum two levels systems) and a classical part is a source of electromagnetic fields 
determining evolution of a quantum part. The general form of a state of quantum part will be 

2"-! 

X = ^i^i where the basic states of it cq, . . . , e2'^-i are simply strings of ones and zeroes of 
the length u where u > n is the size of quantum part which can contain some auxiliary qubits 
behind input for [/, J2 |Ajp = 1, A^ = 2" is the number of all classical input words for U. 

i=0 
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Classical part determines when to "switch on" a tested device (usually it happens many 
times) and when to observe a result of computation. Observation of a state x gives every basic 
state Cj with the corresponding probability |Ajp. 

The problem of recognition of electronic devices presumes the so-called "plug and play" 
technology where a tested device is applied only as a black box. If query transformations are 
only U then our model evidently satisfies requirements of "plug and play" technology where 
we classically control when to switch on a tested device. An implementation of Ucond in the 
framework of this technology is not so easy because it requires a quantum control on applications 
of the device 0. Nevertheless it is possible to implement Ucond in the framework of "plug and 
play" technology. This possibility will be substantiated in the following papers. Now we shall 
simply presume that it is possible. Such difficulty does not exist for the problems of eigenvalue 
and molecular structures recognition. Here we can manage without oracles at all because 
having an explicit form of a quantum gate array realizing a universal function of application 
App we can quantumly control its actions in each element separately and simultaneously thus 
implementing Ucond- 

Let every basic state be partitioned as: 
Cj = I place for code [C], Ri, R2, ■ ■ ■ , Rj), where each register R^ in its turn is partitioned into a 
place for argument, places for time instants and places for the corresponding frequencies. Here 
a complex index i contains one or two integers so that the length of Cj is polynomial of c and 
n of at most second degree. 



2 Obtaining new algorithms from basic quantum tricks 
2.1 GSA and amplitude amplification 

GSA proposed in ( ||Gr|| ) is one of two basic quantum tricks. It is meant for quick getting of 
a quantum state a given the inversion along this state Jq. An inversion along some state d is 
defined by 

T \-\ _ j 1^); if X-La, 
W - j _|^^^ if x = a. 

We also assume that acts like identity if d does not exist. A typical situation is when a 
state is unknown but the inversion along it can be fulfilled easily. Say, let a be a solution of 
equation f{x) = 1 with a simply computable Boolean function /. Then the inversion can 
be implemented by addition modulo 2 of /(x) to an ancillary qubit initialized by This 

transformation maps a state |x, ) to the same state with the sign + or — subject to the 
satisfaction of equality f{x) = 1. The transformation is unitary and can be easily fulfilled given 
a device fulfilling /. All sequential transformations in our formulas will be applied from right 
to left. 

GSA is sequential applications of the transformation G = Ja/g to a state which is chosen 
randomly beforehand. If we apply this transformation 0{^/N) times where is the dimension 

^This would be evidently possible provided we have access to the internal details of our device and can 
quantumly control their work simultaneously. But this assumption contradicts to "plug and play" technology. 
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of main space then an observation of quantum part yields a with visible probability whereas 
without quantum computer we would be compelled to spend of order steps to find a. 

A little difficulty here is that we don't know exactly a time instant t when to stop iterations 
to make probability of error negligible that is needed when applying GSA as subroutine. Here 
the following simple trick helps. 

Define a number B = B{N) so that 1/B is an average value of |(a | 0)| for uniformly 
distributed on a sphere of radius 1 in the space of inputs. A straightforward calculation shows 
that B = 0{\/N). Let GenArg^ be operators generating arbitrary vectors from the space of 
inputs belonging to independent uniform distributions, j G {1, 2, . . . , A;}, and let GenTimeArg^- 
be operators generating time instants tj from independent uniform distributions on integers 
from the segment [0, i?]. Arrange k copies of two working registers: for input and for storage of 
a time instant and fulfill the corresponding operator (jQJa)*^GenArgjGenTimej on each register. 
Now if a exists then the probability to obtain a observing any one register is at least 1/4 (it is 
shown in ||BBHT|| ) and the probability to obtain any fixed other state will be negligible because 
our operators GenArg^ generate independent uniformly distributed samples. If a does not exist 
which means that is identity then the probability to obtain any fixed state will be negligible. 
Denote by the contents of j-th register for argument in the resulting state. Consider the 
following criterion: if at least one fifth of a^, j = 1,2, . . . , k coincide then we decide that a is 
this value, if not then a does not exist. Calculate the error probability of this criterion. Let K 
be a number of such j that dj = d. By the central limit theorem the probability that a fraction 

— (^/'^)~-^ belongs to the segment [ai, 0:2] is closed to -i= / e~^dx. Then the straightforward 



00 ^2 

calculations give that the probability of that K < k/5 will be of order / e~~dx for ai of 

ai 

order ^/k. Thus to make error probability of order 1 / it would suffice to choose k of order 

n = logN. This method can be used not only for GSA but also for other algorithms. If a 

probability to obtain a right result for each of k registers is some positive p which does not 

depend on the dimensionality then to make this probability 1/Ni it is sufficient to choose k 

of order log A^i. In what follows we shall use this simple trick without special mentioning and 

will mark the simultaneous operations of the same kind on all working registers by (g). We 

j 

assume that all ensembles generated by the different j-th copies of operators are taken from 
the independent distributions. 

We shall use the standard norm on operators in Hilbert space defined by \\A\\ = sup ||Aa;||. 

||x||=l 

Given an operator A denote by such operator A' that ||y4 — y4'|| < e. In what follows we 
shall use a simple trick described above making necessary copies of registers and so will raise 
an accuracy of our operators up to required level. When we must repeat an operator T times 
the required accuracy of one application must be 1/T and it may be ensured by only linear 
price in memory as it was shown above. This means that we shall always use instead of A 
without special mentioning where e = 0(1/T) if an operator A must be repeated T times. 



7 



2.2 Revealing of eigenvalues 

The second basic quantum trick is designed to reveal eigenvalues of a given unitary operator 



U. We shall define an operator revealing frequencies in accordance to the work [|AL . 

Let M = 2"*, L = TP . We are going to determine frequencies of unitary operators in within 
1/L where L is a number of applications of V required for the revelation of frequencies with 
this accuracy which means that 1/M is an accuracy that is sufficient to tell apart eigenvalues 
of V . For the recognition of eigenvalues we put p = m + 4 so that L = 16M. 

Denote by (Oi)p a number from [0, 1) of the form l/L. Let an operator U have eigenvalues 
^2mujt: -^j^Qj-Q frequencies u!o,uJi, . . . ,uj]sf'-i are the different real numbers from [0, 1). Denote by 
Ek the space spanned by all eigenvectors corresponding to Uk- An arbitrary vector with the 

N'-l 

length 1 from Ek will be denoted by Thus every state ^ has the form C, = J2 ^k^k- 

k=0 

Let fl = {oJk,i} be some set of integers from {0, 1, . . . , L — 1}, < i < A^ — 1, < k < N' — 1; 
e,6 > 0. Denote L'^{^) = {i : \{O.Uk^i)p — ujk\ < s or 

\{0.ojk,i)p - ^^k - M < £}■ 

Definition 1 A transformation W of the form 

N'-l L-1 
k=0 i=0 

is called a transformation of Ws^s typ^ if for all k and ^ J2 l-^j.fcP > — 25) . 

Thus, 6 is an error probability of getting right frequencies Uk by observation of the second 
register, and e is a precision of frequencies approximations. 

Definition 2 A unitary operator R is called revealing frequencies of U if R belongs to the 
type Wj_^K for any G {1, 2, . . . , L}.| 

The key here is the quantum version of Fourier transform defined by 

. . 1 X T —27Tisl , 

QFT^: \s)-. — J2e^\l) 

1=0 

We need also the following generalization Useq of operator Ucond- 

Ufegk,a) = |U'^x,a). 

This is the result of a sequential applications of U to the main register. To implement this 
operator by means of Ucond fulfill the following cycle. For integer counter j altering from 1 to 
the maximal value L — 1 of a apply f/ iff j < a. Then one cycle consists of Ucond with a properly 
prepared controller and the resulting operator will be U, 
Define an operator revealing frequencies by 

Rev = QFT^ Uf,, QFT^, 



L 

seq' 



^in what follows we shall use this notion only with K = 16. 
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where quantum Fourier transforms are applied to the second register Q In the work [Qz| it 



was proved that Rev is a transformation reveahng frequencies. Now we need more. For the 
redistribution of amphtudes Xk we shall need also a transformation Rest cleaning the second 
register. An ideal candidate for this role would be Rev~^ but a problem is that it requires 
an application of that is physically unrealizable given only device fulfilling U excluding 
evident cases where say U = U~^. We can use this easiest definition of Rest only in case when 
we are given a circuit implementing U (say, gate array) because then is accessible for us as 
well as U. But if C is given only as a black box then the restoring operator should be defined 
separately. 

We shall find an operator restoring ancilla in the form 

Rest = RevD 

where D is some operator of turning. 

Let we are given some integers of the form j;, q - integer, uj^ ^ Uk- Then we could define 
an operator of turning D by D\^k,l) = e"^'^**^-^^^^'''"'' /) where 5k,i = cD^ — (0./)^. It was 
proved in ^ that || (RestRev|x, 0) - |x,0)|| < 7M/L which means that so defined restoring 
operator really restores zeroes in the second register after action of Rev provided L is large 
enough. To create these good approximations we apply a bit more general construction. Put 

D = Enh D Enh 

where an operator Enh calculates an integer function h{l) giving a good approximation {0.h{l))p 
of frequencies in within 1/L given their rough approximation {0.l)m in within 1/M and places 
them into ancilla, D turns each eigenvector on appropriate angle: 

D\^k) = e~^''*(*^~^)((°-''(')^f~(°-'^p^|$fc and the last apphcation of Enh cleans ancilla. An operator 
Enh is accessible given good approximations of eigenvalues. Thus our operator Rest restores 
zeroes in ancilla in within 1/L. 

We can reach an accuracy 1/L of all operators of type Rest that will be less than 1/t 
where t is the number of all steps in computation and this accuracy can be guaranteed with 
logL = p registers. Emphasize that this difficulty with eigenvalue precision arises only when 
U^^ is inaccessible - in the problem of recognition of electronic circuits in the section 3.4 where 
we must choose L ^ M. 

Operators Rev and Rest can be built in the form of quantum gate array using universal 
quantum Klini function App where a code [C] of circuit generating U is a part of input. We 
shall write the operator U corresponding to these two operators as its upper index. 

3 Problems of recognition 

3.1 Getting eigenvectors and recognition of eigenvalues 

Our assumption about sparse spectrum now is stated as L = 16M = const. In view of that 
Rev reveals frequencies it belongs to the type Wj_^^. By the definition of Ws^e it means that 



''the first QFT can be replaced by Walsh-Hadamard transform as in |AL| because on zero ancilla it is 
equivalent 
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Af'-1A/-1 

Rev gives a state J2 J2 K,k\^k, ^k,i) where the seven eighth of probabihty is concentrated on 

fc=0 j=0 

such pairs i, k that {Q.uOi^k)m is closed to uok- This means that we can get eigenvalues with high 
probability observing the second register whereas the first register will contain the corresponding 
eigenvector. This way of getting eigenvectors was proposed in the papers ||AL| , [TM|| . The 
first disadvantage of it is irreversibility. Observing a state we lose full information about it 
and cannot use this state again that is very important for building of nontrivial quantum 
algorithms. The second disadvantage is that this way gives a random eigenvector when it is 
typically required to obtain the eigenvector corresponding to a given frequency. 

Let we are given a good approximation uj^ of some frequency uj written as a string of p its 
sequential binary figures and \et E^j = , . . . $f } be a basis of the subspace E^^ of eigenvectors 

corresponding to all frequencies uj' uj. We shall build an operator State^; which concentrates 

I 

the bulk of amplitude on some superposition of the corresponding eigenvectors: Yl, ^j^'j ^ E^. 

For this aim we are going to apply GSA. Let |a) be some randomly chosen vector from the 
I 

main space: \d) = J2 f^j^^ + J2 i^s^s where all eigenvectors from the second sum correspond to 
j=i 

I 

frequencies uj' ^ uj. Here our target state will be the following vector E^ia) = J2 ^j^^ where 
Xj = —j=^ , that is a vector of the length 1 directed along the projection of a to subspace 

Let A be some set of vectors. We denote by I a an operator changing the sign of all vectors 
from A and remaining unchanged all vectors orthogonal to A. We need to obtain the operator 
Ie^ constrained to the two dimensional subspace S{a,uj) spanned by vectors \a) and E^{a). 

Let ReVj, Restj be j-th copies of the operators Rev, Rest acting on the corresponding places 
of j-th register. Denote by Ij a string contained in the place for frequency of j-th register. Put 

V V 

Ie^ = f^estj Sign^ (g) Revj 

j j 

then Sign^^ changes a sign if and only if for at least 1/2 of all j \{0.lj)p — {0.uj^)p\ < 1/L.f\ 
Applying reasoning from the end of section 2.1 we conclude that the actions of Ie^ and Ie^ 
restricted on S{a, uj) will differ on less than and thus this different can be done very small 
by only linear growth of memory. We thus will omit ~ in our notations. 
We define 

St = GenArg"-^GenTimeArg"^(/5/s^)*GenTimeArg GenArg 

where GenArg and GenTimeArg generate a pair a, [C] and a time instant t correspondingly 
where C is a gate array implementing Jq. Here the actions of la are implemented by the 

^We could choose any fixed p : f < P < | instead of 1/2. Really, so defined Ie^ will change the sign of all 
a € Euj ■ If aJ-Eu; then the probability to obtain ui observing the frequency from Rev is less than i . 
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universal function of application App. Then the result ^ = St|0) of its action on will be 



closed to E^{a). Really, \ {E^{a)\^) \ = | sin(2t arcsin(a|£'j^(a)))| (look at [PBHT|] ). An average 
value of \{a\Ei^{d))\ with the uniformly distributed probability of choice d and t over all space 
and time frame [0, B] correspondingly will be of order 1/ \/N . Thus if t is chosen randomly from 
the uniform distribution over 1, 2, . . . , i? then |(ii^tj(a)|(^)p will have average value not less than 
1/4. Of course it would be much more convenient to obtain E^[d) with the error probability 
converging to zero that is possible by the method described in section 2.1. Namely, arrange h 
equal registers for the states Xfc? k = 1, 2, . . . , in the main space, the corresponding h registers 
for frequencies and associate the variable with each fc-th register. Let St^ be a pattern of St 
operator acting on /c-th register. Remind that operators GenArg;. and GenTimeArg^ generate 
independent distributions for different k = 1,2, ... ,h. Now we define 

State" = Sti (g) St2 (g) . . . (g) St/,. (2) 

So defined operator being applied to zero initial state gives a state Xi X2 <8> • • • <8> X/i where 
an average value of \{E^j{xk)\Xk)\'^ will be closed to some number not less than 1/4 with the 
vanishing probability of error. By the way it means that if we then apply to this state the 
corresponding operators revealing frequencies: Revi (g) Rev2 ... Rev/j then in the resulting 
state X the majority of amplitude will be concentrated on such basic states for which at least 
^ of all registers for frequencies contain numbers I for which |(0./)m — (0.a)^)p| < 1/L. 0. On 
the other hand, if u is not a frequency at all then the probability to obtain such basic state will 
be vanishing in view of independence of distributions generated by GenTimeArg^j, and GenArg^ 
for the different k. 

The time complexity of this algorithm is of order M^/N ■ r? . The last multiplier arises 
due to the copying of registers. Thus we have a solution of the first problem of recognition of 
eigenvalues. 



3.2 Finding thermodynamic functions 

Let we are given a structure of molecule of a gas. The problem is to find its thermodynamic 
function (^. In view of that a summand in this sum quickly converges to zero it is sufficient 
to find few first summands. Thus it is sufficient to be able to find degree of degeneracy of the 
subspace corresponding to frequencies uj' ^ oj for any uj = l/M. Let Eq < Ei < . . . < Eg he 
energy levels of a molecule or eigenvalues of its Hamiltonian H. Then the operator of evolution 
in time frame t will be f/ = e~V*. An addition of a diagonal matrix r ■ I with constant r to 
a Hamiltonian does not change the physical picture. Then choosing r = —Es, t = we 
obtain a unitary operator U which frequencies belong to the segment [0, 1) and thus we reduce 
the problem to the case we have already dealt with. 

Assume that M is fixed and we must examine only few frequencies closed to 0. At first we 
can recognize all numbers of the form l/M that are frequencies in within 1/L. Let u be such 
number. Show how we can find the degree of degeneracy d of the corresponding subspace. This 
is the dimension of subspace E^ spanned by eigenvectors corresponding to frequencies cu' ~ cu. 

^Note that in this criterion ^ could be replaced by any p: 0<p<j-^ — 
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Our plan is the following. Build an operator Ie^ of reflection along this subspace. Then using 
a counting procedure built in the paper PBHT]] we evaluate the time required for turning of 

arbitrary initial vector till this subspace. This time will be about and thus we find d. Fix 
some e > and show how to obtain the value of d in within ed. 

Let operators GenTimeArg^ generate time instants tj from independent uniform distribu- 
tions on the segment [0, [a]], where a is nonnegative number. For a from 1 to \/N fulfil the 
following loop of three steps: 



1. Apply an operator 



{IjEy' GenTimeArg^GenArg 



2. Find the fidelity of result that is the number of all j for which at least ^ — e of all A; are 
such that Uj^k ^ uj. \i the fidelity of this step is larger than on the previous, then proceed 
the loop, if not then stop. 

3. Replace a by 4a/3. 

If by the point 2 we finish computation then the current value a is taken as the rough approx- 
imation of d from above. We have 3a/4 < d < a. To find d more exactly divide the segment 
[3a/4, a] to [1/e] equal parts by points ao < ai < . . . < ai and repeat the procedure from above 
sequentially for all a^. Thus we shall determine d in within g{t)d where function g quickly con- 
verges to zero with e. Thus our algorithm finds d and thermodynamic functions with arbitrary 
relative error in the time 0{\fN)M where the constant depends on admissible error. The more 
refined algorithm can be obtained if we apply the method of counting from the work ||BH'l]] . In 
that work quantum Fourier transform is used analogously to Abrams and Lloyd operator Rev 
only in order to find a time period of function = \G*^^t) that is about JN/d. Their 



method gives the accuracy of order \/d which means that the relative error converges to zero 
iid- 



oo. 



3.3 Recognition of molecular structures 

Now take up a problem of recognition of molecular structures. Here we are given a spectrum 
of molecule and a problem is to recognize its construction. Note that now we have not access 
to a device but it is sufficient to find an arbitrary device generating this spectrum. Clarify 
the formulation assuming the following form of determining spectrum. Let we are given a set 
w = {wi, . . . , Wq} of numbers from [0, 1) of the form Wi = each where G {0, 1, . . . , M — 1}. 
Denote by F a subspace spanned by vectors of the form i = 1, . . .Q. A spectrum S is 
determined by this set w if 

a) for each u & S there exists its good approximation Wi E w: \wi — uj\ < and 

b) each G w is a good approximation of some uj E S. 
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We would obtain a little different formulations of the problem if we want to find a circuit whose 
spectrum only contains one given set of frequencies and/or does not contain some other set, or 
permit some more general form of sparse set for w instead of j^. These versions of the problem 
at hand have similar solutions. 

As above we shall find a recognizing algorithm in GSA form 

(-^O-^cir, w) (3) 

where is arbitrarily chosen vector from the space spanned by codes of circuits, t = 0{VT) 
where T is a number of all possible circuits, and Icir,w is reflection along all such codes [C] that 
Spectr(?7c') is determined by w. Now it is sufficient to build Idr, «>• 

Choose Bf = 0{^/Q) so that a randomly chosen vector w E F satisfies > l/Bf 

with probability 0.99. Let GenFreq^-, GenTimeFreq^- be operators generating correspondingly: 
a linear combination of frequencies ujj G F and time instant t/reg, j < 5/ - all these objects 
from the corresponding uniform distributions over all possible values, and a code of gate array 
generating inversion along the corresponding state Uj. These operators will generate the objects 
in the corresponding ancillary registers. Denote by uj a frequency contained in the j-th register 
(that is initially ojj). 

Assume that a code of circuit generating U is fixed. Define an operator I^r, w by 

Icir, w = ^ GenFreqJ^GenTimeFreqJ^(/BadFreg, wjIcd.Y^'''"'- ' SignGoodFreq 

j 

{hjiBadFreq, w,jf ^ GcuFreq^- GeuTimcFreq^- 

3 

where IsadFreq, uj,j will iuvcrt a sign of states with "bad frequencies" in j-th register that are 
such values of ojj of the form / G {0,1,...,M — 1} which either belong to w and are not a 
good approximation of frequencies uj G Spectr(V") or do not belong to w but have a closed uj G 
Spectr(l^): \ujj — u;| < for all other frequencies this operator acts like identity. Application 
of the sequence preceding SignGoodFreq concentrates amplitude on "bad frequencies". Note 
that Icj can be implemented by a given code by means of quantum Klini operator App. The 
following application of SignGoodFreq inverts a sign of state subject to are there bad frequencies 
or not. Namely, for codes [C] without bad frequencies SignGoodFreq changes the sign, for codes 
[C] with bad frequencies it makes nothing. The following operators clean all ancilla. Thus so 
defined Idr, w will invert a sign of exactly those codes C for which Spectr(f/c') is determined by 
w. We need to define two types of operators: SignGoodFreq and IsadFreq, w,j- 

Associate with each ujj contained in j-th register the family of registers enumerated by two 
indices j, k and containing frequencies Uj^k- 

Definition 3. Call a family of all ujj^k good if for at least 1/5 from all j the following 
property takes place: for at least 1/10 of all k ujj^k ^ ujj E w. 

Registers enumerated by different k for a fixed j are designed for the application of j-th 
copy of operator State'^ defined in the previous section. Here it has the form State'^-'. Each k 
corresponds to the operator St^ from the definition (Q) so that each uj^k will be a frequency 
obtained from the result of St^ action. 
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At first build an operator lBadFreq,w,j- Put 



BadFreq,w,j 



(State'^0 ^Restj-fcj Sign' (g) [Rev j-^State'^ 



j,k j,k 



where an operator Sign' will change a sign of only states with bad families of frequencies. 

If a frequency Uj is bad then in the previous section it was shown that only for vanishing part 
of all k we can have uj^k ^ ujj E w and before Sign' almost all probability will be concentrated 
on bad families Uj^k, hence lBadFreq,w,j changes the sign. 

If ujj is good then it belongs to w and has a closed u' G S. By the previous section for 
about I ■ ;| = ^ > I of all will be Uj^k ~ u; e w and before Sign' almost all probability will 
be concentrated on good families, hence the sign will not be changed. 

Thus lBadFreq,w,j is defined correctly. 

Put 

SignGoodFreq = (g) [(State'^O^^Restj- J Sign(g) [Rev^- ^State*^^] 

j,k j,k 

where an operator Sign changes a sign only for states with good families of frequencies. If a 
frequency ujj is not bad then for about ^ • | = ^ of all will be ujj^k ~ E w. If a frequency 
ujj is bad then we can obtain ujj^k ^ ujj E w only for the vanishing part of k as it was shown in 
the previous section. Thus SignGoodFreq acts how it is needed[} 

Now calculate the complexity of our algorithm recognizing molecular circuit. The first 
multiplier \/T issues immediately from (j^). The following multiplier as y/Q issues from the 
immediate definition of Idr, w- At last the definition of iBadFreq, w brings the multipher M^/N. 
The resulting complexity will be of order My/TNQn"^. 



3.4 Distinguishing of eigenvectors of two operators with the same 
eigenvalue 

Now we are going to take up the most difficult of our problems - a problem of recognition of 
electronic devices. The difficulty is that here we need not to find a circuit with given spectrum 
but to simulate an action of given circuit. Remind that now we assume that frequencies can 
be determined in within 1/L given their approximation in within 1/M where L 3> M. 

As a first step we consider the following question. Given two operators U and V having 
the same eigenvalue u how to find a difference between the corresponding eigenvectors? Let 
L^, be the subspaces spanned by eigenvectors of U and V corresponding to all frequencies 
cu' ~ tu. (A particular case is when u; is a frequency of U but not of V . Here L^J^ = and 
our algorithm will work in this situation.) We shall omit the index uj in the notations. For 
u G L^, = 1 put 



Hu = minlyl — |('u|f)P | v G L^, ||f || = 1}. 

^ Again we could take arbitrary pi : < pi < 1 instead of and p2 '. < p2 < instead of 1/5 in the 
definition of a good family. 
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that is the sine of angle between a vector u and the subspace or a distance between u and 
this subspace, and analogously define /ly for v e , \\v\\ = 1. 
Put 

/lu — max //y = max /ly. 

Then say iijj = Q means that U QV . li the dimension of spaces LF , LX are equal then /lu = /ly, 
if they are not equal, say dimL^ > dimL^ then ixjj = 1. 
Let d — d{N) be some function taking values from (0, 1]. 

Call these subspaccs d- distinguishable if some of /xy is not less than d, or some of the 
subspaces is empty and another is nonempty. 

We shall build a procedure that determines are these subspaces the same or not provided 
they can be either d- distinguishable or coincident. The less the function d{N) is the more 
accurate our recognition will be. Let LF nlX — Lq. Then LF — Lq® L'jj and = Lq L'y. 
Note that if L'^ ^ then for all vectors from of the length 1 their distances from are 
exactly and the same thing takes place with if L'y is not empty. Let L' be linear subspace 
spanned by vectors from L'yVJL'jj. Denote by Proj aB a projection of subspace B to subspace A. 
If dim > dim then we have the following expansion to the sum of orthogonal subspaces: 
L'^ — L'/r0Proj luLF 1 where L'(j is a subspace in consisting of vectors orthogonal to LF . 
Let Ly be defined symmetrically. 

Then either 

. = or 

• dim = dim and L' ^ 0, or 

• dim > dim and L'lj 0, or 

• dim < dim and L'{, ^ 0. 

We define the main operator determining the equality of iF and LF by 

Difference = Differ"^ SignDif Differ, , , 

Differ = Dif^^^e dimDifif >LvDifLt^<L^Difrt^>LvDif^u<iV 

where SignDif changes a sign of main ancilla a^ij iff at least one ancilla from the list a — 
{(^same dim,(^L'^>L'^,OiLU^LV,a'^uyLV,o:'^u^Lv} contains 1, and each operator of the sort Dif 
changes the corresponding ancilla from a in cases 

• dim = dimL^ and iF ^ , 

• dim LF > dim and nv < ^2/3, 

• dim < dim and /lu < -/s/S, 

• dimL^ > dimL^ and nv > \Jl/?>, or = 0, 
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• dim < dim and /lu > 073, or = 

correspondingly and all of them do nothing if — . In view of symmetry it is sufficient 
to define Dif operators in the first, second and fourth cases. Note that the first case dimL*^ = 
dimL^ is the only non-degenerate case and the definition of Dif here will be the most difficult. 

Definition of Dif^a^e dim- 

Suppose that dimL'^ = dimL^. Our first aim now is to build an operator Inv which acts 
as identity if and are coincident and acts like I^/ if they are d- distinguishable. Arrange 
the first two ancillary qubits au, O-v which will signal that a state at hand has a projection at 
least of the length 1/3 to or, correspondingly to . Consider the following operator 



Check = Rest^Ancy (g) Rev^ (g) Rest^ Anc[/ (g) Rev^ 



where Anc inverts the corresponding ancilla if and only if at least nine tenth of copies for the 
respective frequencies are equal to uj in within 1/M. It coincides with the inverse operator 
Check-\ 

Let t be some random integer from the segment [0, | ]. We define the following operator 
of Grover's type: 

Turn* = {hulLvf (5) 

Call two subspaces and LX almost orthogonal iff for some G \/l — ji^ < 1/30. 

If L'^ and are not almost orthogonal then given some a e L^ (a G Ly) an average distance 
between Turntja) and Lu (Ly) will be at least 1/2 if and are d- distinguishable and zero 
if these subspaces are coincident. To tell apart close location and almost orthogonality build 
two operators: Distort and DistcZosed- 

At first suppose that and are almost orthogonal. Then au = I means that ay = 0. 
Introduce a notation 

r L^, if au = 1, 
1 L^, if av = 1. 



L{au, ay] 



Let a be a vector at hand from the space of inputs. Note that if ^ then for each a_LL' 
we shall have au = ay because such a belongs to the subspace spanned by Lq and orthogonal 
subspace to LF U . The first operator Distort will do nothing if au = olv and will change a 
sign and a special ancilla aort if a projection of a to L{au, cty) is less than 1/30. 

The second operator Distc^osed will act like identity if au — OLy and will change a sign 
in case when the following conditions are satisfied simultaneously: a & L' , LF and LX are 
distinguishable, a^rt — 0. 

Put 

Distort = (g) RcSjSi^^ (g) Rcj 

where Re (Res) denotes Rev^ (Rest^) if au — l,ay — 0, Rev*^ (Rest*^) if ay — l,au — 0, 
and identity if au — (Xy] Si^^ changes a sign inverting simultaneously aort iff at least a half 
of frequencies cuj are such that \u!j — lu] > 1/M and au Oiy- If we want to clean the second 
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ancilla after the action of Distort and remain change in sign then we can use an operator 
Dist~^j = ReSjS^^^ {g)^ Re^ where S acts hke Si only without changing a sign. 
The second operator will be defined by the following equations 



DisUosed = D^' . . . D-'S'Dr,Dn-i . . . Di, 

Dj = (GenTimeArgj.)-i(Turn(.)-^ ORestJ^^ 
j = l,2,...,n, 



(8)Rev^, 



3,' 



Turn^ . GenTimeArg ■ 



where operator Sig^^ changes the corresponding ancilla Pj only in one of the two cases: 

1) Of/ = 1 and at least a half of cuj^k are such that \u)j^k — (^1 > '^/M, or 

2) au = 0, av = I and at least a half of cuj^k are such that \uj^k — uj\ < 1/M. 

An operator 5"' changes a sign iff some of au, ay is nonzero and at least 1/20 of all Pj contain 
1. 

Consider the action of Distdosed following to Check on an input vector a. Let at first 
L'^ ^ which means that they are distinguishable. 

If a_L , LX then au — OLy = ^ and Distdosed makes nothing. 

If a G Lq then au = ay = 1 and all Sig-^^ makes nothing because for almost all j about 
3/4 of ujj^k are closed to uj: \u!j^k — ^| < '^/M, hence S' and Distdosed do nothing. 

Let a e L'. Prove that Distdosed changes a sign. Expand L' to the sum of orthogonal 
subspaces: L' — L'jj^L'^'^*. Denote the result of action of Turn^^. on a by 

If q; G L'jj then a^/ = 1, and for more than 1/10 of all revealed frequencies are not closed 
to oj with probability about | ■ hence a sign will be changed by the point 1). 

If a G L'u°^^ then by the same reason we obtain the change of sign by the point 2). Hence 
Distciosed changes a sign for all a G L' . 

Now we can define Inv: 

Inv = Check Dist^^^Distc/osedDistort Check. 

For a_L LF , we have Inv|a) = \a) because Check gives zero in ancilla au-, ay thus depriv- 
ing the following operators ability to change somehow a state vector. If a G Lq then Inv|a) = \a) 
because Distort makes nothing and Distc^osed makes nothing as well. Thus Inv|a) = \a) for a±L', 
and Inv|a) = —\a) for a G L'. 

Now we arc ready to build an operator Dif^^^e dim inverting the ancilla agame dim if and only 
if L'^ and are distinguishable. Let Gen generate a list y, [ly], [Cz] where [Cz] is a code of 
circuit generating some unitary operator Z = having only eigenvalues 1 and —1 that is its 
frequencies are and 1/2 and the space corresponding to frequency is one dimensional where 
y is its basic vector. As usually index j means that the corresponding vectors yj are taken from 
the uniform distribution on all possible vectors. Assume that operators of the form Gen~^ are 
accessible for us as well. Put 

Difsome dim = <S> GenTimeArgJ^Gen~^(InVj ly.y^Restf' Change 

z (6) 
ReVj^{Iy. Invj)*j GenjGenTimeArgj 

j 
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where each copy of Inv acts on the register where initially is placed yj , Change makes a desired 
change in a resulting qubit asame dim provided at least 5/32 of all frequencies differ from in 
more than 1/M. 

The group {ly. InVjY^ of GSA type turns essentially a vector yj generated by Gen^ if and 
only if and are d- distinguishable. 

If L'^ = LX then yj remains unchanged and at least 7/8 of all frequencies will be closed to 

0. 

If LF 7^ LX then for the result of the turn of yj at least | ■ ^ = ^ of frequencies will be far 
from because they must be closed to 1/2. 

Definition oi'Diiiu^iv . 

Suppose that dimL^ > dimL^ and /ly < y^2/3. Remind that here we have an expansion 
to the sum of orthogonal subspaces LF = L'lj 0Proj lvLF where L'^ ^ 0. We shall define the 
operator Dif by a very similar way as in previous case: 



Dific/>iv = (g) [GenTimeArgTiGenTi(Inv",, /„J*^Restf^ 



RevJ^(/y^ InvJ;y)*-'GenjGenTimeArg^- 







Change 



where the definition of Inv'^ inverting L'^ looks like Distort only L'^ will play a role of V: 



Inv'^ 



Check 









(g)Resr 


Si^w 


(g)Rer 


. k 




. k 



Check. 



^~ V 



~ V 



Here Re' and Res' act like Rev^ and Rest^ only ii ajj = 1 and if «[/ = then they do 
nothing, Si^^ changes a sign only in one case: if a^/ = 1 and at least 3/4 of all frequencies uj^ 
are far from uj: Icj, — c^l > 1/M. Thus in Dif operator we shall use a set of ancillary registers 
enumerated by pairs of indices j, k. 

For (ij G Proj luLX in view of /iy < y^2/3 an operator Si^^^ does not change a sign because 
here the fraction of all frequencies closed toa;is|-| = ^>|. 
For Oj-LProj l^L^ an operator Inv'^ makes nothing. 

Definition of Bif^S^^v 

Suppose that dimL^ > dimL^ and nv > v 1/3. The definition of Dif will be similar to the 



previous case only the whole subspace Lu will play a role of L': 

Dir^i^LV = ® [GenTimeArg-^Gen7^(InVj-C7 ^J^^Rest/ 



Change 



Rev^-^(/j^^. InVj_[/)*-'GenjGenTimeArgj 

j 



where 



Inv 



u 



Check 









(g)R~esr 




(g)Rer 


. k 




. k 



Check. 



^Tlius wc could take any number p 



< P < instead of in the definition of Change. 
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Here Si_^^ changes a sign if more than a half of frequencies are far from uj: \ujj —uj\> 1/M. The 
satisfying of the conditions required for Dif operator is based now on inequahty || = ^ > | 
and can be checked straightforwardly. 

At last estimate the complexity of constructed procedure. An operator Turn requires 
of order Tunicompiexity = M^l/d elementary steps. Then, Difference (|) requires of order 

Tmiacompiexity''/N that is 0{M^ N/d) elementary steps. Note that there exists the similar form 
of operator Difference which does not act on resulting qubit adif but changes a sign instead and 
such operator can be constructed similarly. Denote this operator by Difference^jgn. Assume 
that an input of it contains a frequency uj. 



3.5 Recognition of electronic devices circuits 

Now we are ready to take up the recognition of circuits. We assume that for every pair of 
circuits for their transformations Ui,U2 subspaces spanned by corresponding eigenvalues are 
either coincident or d- distinguishable. Assume also that our coding procedure gives one- 
to-one correspondence between circuits and T basic states eo, ei, . . . , er-i in the space Hdr- 
Recognizing procedure is denoted by Rec and will have the GSA form: 

Rec = (Vt/)*, t = 0{VT) (7) 

acting on states of the form \x) where basic states for x are codes of circuits. Here G H^r is 
chosen arbitrarily and lu inverts a sign of every code which circuit induces a given operator U. 
An implementation of /q is straightforward and all that we need is to build lu- 
We define lu as 

ConcjY^g^^ jDifferencCj Sign(^ [Difference jConc/reg, j] 

j 

where Cone freq for every basic state C of argument will generate some arbitrary distribution 
of amplitude on ancillary register with Q basic states and then will concentrate substantial 
part of amplitude on a frequency uj for which L'^ and are distinguishable (if such frequency 
exists). Then operator Difference^ changes resulting qubit for jth copy if and only if on this 
frequency these subspaces are distinguishable. The following operator Sign changes a sign if 
and only if at least one fifth of resulting qubits adif contain 1, e.g. if and only if operators U 
and Uc are the same. Then the following applications of DifferencCj to each copy of register 
clean the corresponding resulting qubits and inverses operators for Concj restore an initial state 
of ancillary register. Difference was constructed in the previous section and all we need now is 
to build Cone freq, j- This transformation can be defined as 



Cone freq, j = GcuTimeFreq^ ^ GenFreq^ ^ {DiSerencesignlivjY^ GenFreq^ GenTimeFreq^- (8) 

If U and Uc are different then by our assumption for some u their subspaces and 
are ^-distinguishable then Concj will concentrate substantially large part of amplitude over all 
j on some combination of such values u. Thus we have constructed a required procedure Rec 
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which gives a target code with substantial probabihty as a result of observation of the register 
for code C . After the observation we can verify a fitness of a code found C by a straightforward 
procedure. It is similar to lu with the single change: Sign will be replaced by change in a 
special ancilla which can be observed after procedure and thus we shall learn does a code C at 
hand fit or not. 

What is the complexity of our procedure Rec? The complexity Mn^ ^ N / d of Difference 
must be multiplied by y/Q issued from (|^) and by \/T issued from the immediate definition 
(0). The resulting complexity will be Mn^ JTQN / d. 



3.6 Advantages of the recognizing algorithms 

Advantages of the proposed algorithms are connected with their high speed and small memory. 
Particularly, the algorithm for molecular structures recognition makes possible to recognize 
molecular circuits using microscopic memory whereas classically this task requires exponentially 
large memory. Compare the proposed algorithms with their classical counterparts. We shall 
omit logarithmic multipliers. 

1. Recognition of eigenvalues and finding thermodynamic functions. Fix some value of 
M determining a precision of eigenvalue approximation. Consider at first the case when the 
number of ancillary qubits in a quantum gate array at hand is small. Then by the direct 
classical method we must build a matrix of unitary transform induced by a gate array. It 
requires of order A^^ steps and at least of order A^^ bits. The known quantum algorithm 
given by Travaglione and Milburn in |[rM|| basing plainly on Abrams and Lloyd operator Rev 
contains repeated measurements of frequencies hence it requires the time of order NM - for 
sparse spectrums it is of the same order as for Hams-Raedt algorithm and its only advantage 
over the last is exponential memory saving. 

Our algorithm recognizes an eigenvalue in \fNMn steps. This time for the sparse area of 
spectrum is about square root of the time of best known algorithms. Here the memory will 



be of order g"^ qubits {g is the size of gate array), that is about squared memory used in ||AL 
but still exponentially smaller than of classical methods. Thus the proposed algorithm gives 
essential speedup over known methods in case when the number of ancillary qubits in a given 
gate array is small (as in case of molecular structure simulated by gate array) and an area of 
spectrum at hand is sparse. The same advantage we have with the proposed method of finding 
thermodynamic functions. 

If spectrums are dense we assume that M = N which means that eigenvalues differ at least 
on 1/A^. Then the time of our algorithm is 0{N). 

Consider the case when the number a of ancillary qubits involved simultaneously in the gate 
array is much more than the length n of input. Then the direct classical method requires more 
than 2^" steps and at least 2"^ bits whereas our algorithm requires only about g2'^ steps and 
gn"^ memory and the quantum speedup may be more than square root. 

2. The recognition of molecular structures. At first assume that spectrums are sparse. To 
be able to compare our method with the evident classical algorithm let us assume that a code of 
molecular circuit of the length n is a string of ones and zeroes of this length. Thus M = N. The 
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next natural assumption that may be also presumed for electronic circuits is that the sampling 
of a code of circuit from the uniform distribution induces a sampling of all possible spectrums 
from the uniform distribution as well. Then the number of all possible choices of spectrums 
approximations (or parts of spectrum subject to the statement of recognition problem) in within 
1/ L consisting of frequencies of the form //M is about 2^ = A^. It means that in our assumption 
M and Q must be logarithmic of A^. Hence our method has the time complexity 0{N). With 
these assumptions the time complexity of the classical direct algorithm examining all codes and 
calculating the corresponding spectrums is about ■ N = whereas our algorithm requires 
the time about and logarithmic memory. Thus the quantum time for this problem is about 
fourth root of the time of classical direct method and quantum space is logarithmic. 

If spectrums are dense then Q and M will be of order and our method requires the time 
0{N'^'^) comparatively with 0{N'^) of direct classical way. 

3. Recognition of electronic devices. Here in the general case there are no classical analogs. 
Compare two algorithm constructed above with their classical and known quantum counter- 
parts. At first consider the single recognizing quantum algorithm that can be easily deduced 
from the technique known before. This is an algorithm of recognizing a circuit realizing clas- 
sical involutive function of the form / : Q — > Q, f = f~^. This task can be reduced to 
the search of y such that the following logic formula is true: Vx A{x,y) where A{x,y) is some 
predicate. Indeed, if we take Y{x) = U{x) in place of A{x, y) where F is a function whose code 
is y then we just obtain the problem of recognition of circuit generating U. An algorithm for 
such formulas given in ||BCW|| has the time complexity of order yTN . This task is a particular 
case of our algorithm for involutive devices and it has the same complexity. In this particular 
case quantum time is of order square root of classical. But if we regard a bit more general 
but still restricted problem of recognition of involutive devices producing linear combinations 
of basic states (like quantum subroutines) an advantage over classical method of recognition 
will be more. For example, consider the restricted problem when we must choose between 
two alternative constructions of a tested device inducing not classical unitary transformation. 
The naive method of observing the results of action of a tested device on the different inputs 
requires of order ^A^^ steps to restore the matrix of the operator Uc in within e. Then this e 
must be less than 1/ \fN to give vanishing difference between operators in Hilbert space. Hence 
the time complexity of the naive method of recognition is roughly A^''/^. On the other hand 
the method proposed in the section 3.4 requires the choice of d only converging to zero with 
A^ converging to infinity. Thus the time required by our method is a little more than -y/iV. We 
thus have almost seventh degree speedup for the problem of distinguishing electronic circuits 
generating transformations with not classical matrices. 



4 Conclusion 

The main conclusion is that molecular structure and physical properties of environment can 
be quickly recognized on the microscopic level whereas the classical methods require huge time 
and especially memory. The new algorithms recognizing eigenvalues with fixed precision and 
molecular structure, finding thermodynamic functions give a quadratic speedup comparatively 
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with the best classical algorithms and exponential memory saving. The new method based on 
quantum computing was proposed for fast recognition of electronic devices. By this method 
two devices with the same given spectrum can be distinguished in the time about seventh root 
of the time of direct measurements. All these algorithms show essential potential advantages of 
microscopic sized quantum devices comparatively with their classical counterparts with much 
bigger memory. The advantages touch intellectual tasks like recognition of the structure of 
other devices and important properties of environment. The proposed algorithms are built of 
standard known subroutines; they have simple structure and lay completely in the framework 
of conventional paradigm of quantum computing. 
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